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(57) Abstract: The invention concerns a seismic prospecting method which 
consists in emitting underground a compression^ seismic wave and picking up 
with sensors seismic data having at least a shearing component and in process- 
ing the data corresponding to said shearing component to deduce therefrom an 
information concerning the underground geological structure. The invention is 
characterised in that the method consists in determining an estimate of the ratio 
(a), wherein vp and vs are real local speed values of compression and shearing, 
wherein 1 represents the underground depth co-ordinate, and wherein z is the 
value of said depth co-ordinate at the bottom surface of the last layer to be anal- 
ysed, and by inverting the seismic data so as to deduce the local speed values 
of compression and shearing for said layer to be analysed, using a modelling 
process wherein said estimate is used for the first invariant parameter y cff . 

(57) Abrege: Procede de prospection sis mi que selon lequel on emet dans le 
sous-sol une onde sismique de compression et on recueille a l'aide de capteurs 
des donnees sismiques ayant au moins une composante de cisailiement et selon 
lequel on traite les donnees correspondant a cette composante de cisailiement 
pour en deduire une information sur la ge*ologie du sous-sol, caracterise en ce 
qu'on determine une estimation du rapport (a), ou vp et vs sont les valeurs de 
vitesses reelles locales de compression et de cisailiement, ou 1 dCsigne la coor- 
donnee de la profondeur dans le sous-sol, et ou z est la valeur de cette coordon- 
nee de profondeur a la surface inferieure de la demiere couche a analyser, et on 
inverse les donmSes sismiques afin de deduire les valeurs de vitesses locales de 
compression et de cisailiement pour ladite couche a analyser, en utilisant une 
mode'lisation dans laqueile cette estimation est utilised pour le parametre inva- 
riant Ycff. 
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PROCEDE DE PROSPECTION SISMIQUE METTANT EN OEUVRE UN 
TRAITEMENT SUR LES ONDES CONVERT1ES 

La pr6sente invention est relative a un proc6de de prospection 
5 sismique rnettant en oeuvre un traitement sur les ondes converties. 

Le principe g6n6ral de la prospection sismique consiste a 
provoquer, d I'aide d une source sismique, un 6branlement dans le sous-sol 
et £ enregistrer, a ('aide de capteurs, des donn6es sismiques generees par 
cet ebranlement pour en tirer une information sur la geologie du sous-sol et 

10 en particulier pour d6tecter la presence d'hydrocarbure. 

On a represents sur la figure 1 une onde acoustique se propageant 
dans un sous-sol a partir d'une source 1 . Cette onde acoustique est, dans 
I'exemple represent6, une onde de compression qui se reflechit dans le 
sous-sol en se d6composant selon une onde reflechie de compression et 

1 5 une onde r&fiechie de cisaillement. 

On rappelle que les ondes de compression (ondes dites de type P) 
vibrent dans leur direction de propagation, tandis que les ondes de 
cisaillement (ondes dites de type S) vibrent quant a elies 
perpendiculairement a cette direction de propagation. La Vitesse de 

20 propagation des ondes de cisaillement est inferieure a la Vitesse de 
propagation des ondes de compression et la connaissance des champs de 
vitesse des ondes de compression et des ondes de cisaillement permet de 
determiner les informations sur le sous-sol. Par exemple, le rapport entre la 
vitesse des ondes de compression et la vitesse des ondes de cisaillement 

25 permet de determiner le coefficient de pression des roches traversees et 
sert egalement d'indicateur de presence d'hydrocarbure. 

Classiquement, pour inverser des donnees sismiques, on utilise des 
modelisations des champs de vitesses qui dependent d'un certain nombre 
de parametres - que Ton considere comme etant invariants pour une 

30 gamme de deports source/recepteur donnee et une zone spatiale 
d'acquisition donn6e t ces parametres pouvant neanmoins etre "lentement" 
variables spatialement, c'est S dire etre variables d'une zone spatiale a une 
autre. 
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Pour inverser des donnees sismiques correspondant a des 
reflexions SS ou PP dans le sous-sol, on utilise des parametres Vp et Vs 
qui represented des vitesses apparentes des ondes de compression et des 
ondes de cisaillement apres correction dynamique ("Normal Move Out" ou 
5 "NMO" selon la terminologie anglo-saxonne generalement utilisee par 
I'homme du metier), ainsi que des parametres Tp et Ts t qui represented 
respectivement des temps de trajets verticaux des ondes P et S. Les 
parametres Tp et Vp suffisent pour I'analyse des vitesses PP, tandis que les 
parametres Ts et Vs suffisent pour Tanalyse des vitesses SS. 
10 L'analyse de vitesses converges (reflexions PS) se fait quant a elle 

g6neralement en utilisant, dans le domaine temporel, des moderations 
faisant intervener les parametres Vp, Vs, ainsi qu'un parametre Vc, ou Vc est 
tel que : 

Tc.Vc 2 = Tp. Vp 2 + Ts. Vs 2 , avec Tc = Ts+Tp. 
15 Les moderations utilisant ces trois parametres sont performantes 

dans le cas de materiaux homogenes et isotropes pour les ondes S et P. 
Par contre dans le cas de milieux verticalement inhomogenes ou a forte 
anisotropie, il a 6te montre qu'il convenait de tenir compte de deux autres 
parametres, d6sign6s dans la litterature par y efT et y 0 . ou y eff = y n 2 /yo avec 
20 y n = VpA/s et yo = Ts/Tp. 

Le deport xc du point de reflexion par rapport a une source depend 
en effet au premier ordre du paramfetre y e ff et au second ordre du parametre 
Yo, ainsi que de la quantite Tc.Vc 2 . 

On pourra a cet egard avantageusement se referer a la publication 
25 suivante : 

- [1] Thomsen, L, 1998, "Converted-Wave reflection seismology over 
anisotropic, inhomogeneous media". 68th annual meeting, SEG Expanded 
Abstracts, 2048-2051. 

Toutefois, dans cette publication, le parametre Yeff est suppose 
30 connu. Or, dans la pratique, aucun des parametres precites n'est 
habituellement imm6diatement connu. 
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Uinvention a quand & elle pour but de proposer un procede de 
traitement sismique sur les ondes converties qui est particulterement fiabie 
et ce ind6pendamment de la connaissance pr6alable des paramfetres y e ff et 
Yo. 

5 R6cemment, il a ete propose de determiner le d6port lateral du point 

de conversion en utilisant la correlation laterale entre d'une part des images 
a d§port source/recepteur avant et d'autre part des images a deport 
source/recepteur arriere, c'est a dire des images obtenues en inversant les 
positions des sources et des recepteurs. 

10 On pourra a cet egard se r6ferer a : 

- [1] Herrmann, P., Michaud. G., Granger P.Y., 1999, "Stacking mode- 
converted waves", presented at the CSEG conference, Calgary, May 1999. 

Uinvention propose quant a elle un procede de prospection 
sismique selon lequel on emet dans le sous-sol une onde sismique de 

15 compression et on recueille £ Paide de capteurs des donnees sismiques 
ayant au moins une composante de cisaillement et selon lequel on traite les 
donn6es correspondant a cette composante de cisaillement pour en deduire 
une information sur la geologie du sous-sol, caracterise en ce qu'on 

determine une estimation du rapport j* vpAl ^/ ^vs.dl, ou vp et vs sont 

20 les valeurs de vitesses r6elles locales de compression et de cisaillement, ou 
I d£signe la coordonnee de la profondeur dans le sous-sol, ou z est la 
valeur de cette coordonnee de profondeur a la surface inf6rieure de la 
derniere couche a analyser et ou z 0 est la valeur de cette coordonnee de 
profondeur £ la surface superieure de cette couche ou d'une couche au 

25 dessus de celle-ci, 

et on inverse les donnees sismiques afin de deduire les valeurs de vitesses 
locales de compression et de cisaillement pour ladite couche a analyser, en 
utilisant une modelisation dans laquelle cette estimation est utilisee pour le 
parametre invariant y e tf. 

30 ^invention est avantageusement cornpletee par les drfferentes 

caracteristiques suivantes prises seules ou selon toutes leurs combinaisons 
techniquement possibles : 
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- on determine le parametre y^fr en mettant en oeuvre, pour differentes 
valeurs possibles pour celui-ci, un traitement de migration des donnees 
sismiques qui correspondent d la composante de cisaillement, et en 
d6terminant la valeur du parametre y eff pour laquelle les images sismiques 

5 avant et arriere sont les mieux corr6lees ; 

- pour faire varier le parametre y efft on pose 

vp a = a vp 0 et vs 3 = p vs 0 , 

ou vp 0 et vs 0 sont les valeurs approximates prealablement 
determines pour vp et vs, 

10 et on fait varier a la fois les variables a et p ; 

- la mod§lisation fait intervenir comme parametres invariants au moins 
quatre des parametres y 0 , Yeff, Tp, Fp, Tc, Fc, avec y 0 = TsHp, y eff = Fp/Fs, 
Tc = Tp + Ts f ou Tp et Ts represented respectivement des temps de 
trajets verticaux des ondes de compression et de cisaillement, ou Fp est tel 

15 que (Fp/Tp) 1/2 represente une Vitesse de compression et ou Fc est tel que 
((Fc-Fp)/Ts) 1/2 represente une Vitesse de cisaillement. 

- lorsque Ton fait varier les variables a et 0, on remplace les parametres y 0 , 
yeff, Tp, Fp, Fc par 

Yo' = a/p * y 0 
20 yeff 1 = a/p * yeff 

Tp' = TpM1+Yo)/(1+yo') 

Fp' = Fp* a 2 M1+Yo)/(1+yo') 

Fc' = Fc * ap * (1+y 0 )/(1+ Yo') * (1+ y* f J/(1+ y«n) 
et on met en ceuvre sur les donnees sismiques une migration 
25 correspondant d ces nouveaux parametres ; 

- pour faire varier le parametre y eff , on pose p = 1/a et on fait varier a ; 

- a Tissue de la determination du parametre yeff, on fait varier vp et vs tout en 
maintenant y e ff constant et on determine le parametre Fc pour lequel 
I'alignement selon la direction de deport est maximal ; 
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5 

- pour faire varier vp et vs, on pose vp a = a vpi et vs a - a vsi , ou vpi et vsi 
sont des valeurs dtterminees pour vp et vs dans Tetape 2, et on fait varier la 
variable a ; 

- a Tissue de la determination du parametre Fc, on determine le paramfetre 
5 Tp et/ou le parametre y 0 = Ts/Tp ; 

- ie parametre Tp est avantageusement determine a parti r du champ de 
vitesse vp determine a partir de I'analyse de la composante de compression 
des donnees sismiques ; 

- on met ensuite en oeuvre un traitement de mise en profondeur commune 
1 0 des modeles de vitesses S et P ; 

- a Tissue du traitement de mise en profondeur, on met en oeuvre un 
traitement a courbure a grand deport en faisant varier les parametres 
d'anisotropie 8 et a tout en maintenant constant le rapport (1+25)/(1+2a). 

' D'autres caracteristiques et avantages de T invention ressortiront 
15 encore de la description qui suit. Cette description est purement illustrative 
et non limitative. Elle doit etre lue en regard des figures annexees sur 
lesquelles : 

- la figure 1, deja discutee, illustre schematiquement la decomposition en 
ondes S et P d'une onde de compression emise par une source 

20 acoustique ; 

- la figure 2 illustre differentes 6tapes de traitement conformes a une mise 
en oeuvre possible pour Tinvention ; 

- les figures 3a et 3b illustrent Tinfluence du positionnement du point de 
conversion commun dans la determination de y e ff. 

25 Le traitement qui est decrit ci-apres en reference aux figures 2 et 

suivantes se met en oeuvre sur des images sismiques constitutes par des 
traces sismiques correspondant les unes a des ondes de type PS, les 
autres a des ondes de type PP, acquises au moyen d'un ensemble dense 
de sources sismiques et de capteurs geophones ou hydrophones 2D ou 3D 

30 repartis selon au moins une direction d'acquisition. 

Ce traitement comporte diff6rentes 6tapes referencees de 1 a 6. 
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Dans une premiere etape (etape 1), on determine des valeurs 
approchees pour les parametres de moderation des vitesses dans la 
couche d'analyse, en principe la premiere couche du sous-sol dont les 
vitesses sont a inverser, les vitesses des couches au dessus de cette 
5 couche d'analyse ayant quant a elles ete prealablement inversees, par 
exemple au moyen du meme traitement. 

Notamment, on determine une estimation du parametre y e ff. 
A cet effet, on utilise par exemple les images sismiques ternporelles 
correspondant d ('acquisition PP initiate pour determiner une valeur 
10 d'initialisation pour le parametre invariant Vp et on utilise les images 
sismiques ternporelles correspondant a I'acquisition PS initiate pour 
determiner une valeur d'initialisgtion pour le parametre invariant Vs. 

Dans un deuxieme temps (etape 2), on met en ceuvre un traitement 
destin6 £ determiner une valeur de y e ff plus precise. 
15 A cet effet, on determine une valeur de y e ff pour laquelle le decalage 

AX, selon une direction spatiale d'acquisition, entre les images sismiques 
PS avant et les images sismiques PS arriere - lesquelles sont obtenues par 
inversion des positions des emetteurs et des sources le long de la direction 
d'acquisition consid6ree - est nul. 
20 Plus pr6cisement t les images sismiques avant sont obtenues en 

consid6rant un ensemble dense d'emetteurs alignes suivant Tazimut qui 
correspond a la direction d'acquisition et en relevant les traces sismiques 
obtenues pour les capteurs d6caies par rapport aux emetteurs, dans la 
direction d'acquisition, d'un deport (algebrique) constant pre-defini. Les 
25 images arriere sont obtenues en utilisant les emetteurs directement voisins 
des capteurs utilises pour Tacquisition des images sismiques avant, ainsi 
que des capteurs decales par rapport a ces emetteurs, dans la direction 
d'acquisition, d'un deport (algebrique) inverse de celui utilise pour 
I'acquisition des images avant. 
30 Ce choix d'emetteurs et capteurs est extrait des donnees lors de la 

migration. 
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Les images avant et arriere ainsi determinees sont ensuite utilisees 
pour determiner la valeur correcte de y cff . 

En effet, on a vu qu'il avait §t§ montre que le deport xc du point de 
conversion par rapport £ la source, dans la direction d'acquisition, depend 
5 de yeff au premier ordre et du parametre 70, ainsi que de la quantity Tc.\/c 2 
au second ordre. 

Or, la determination de cette quantite Tc.Vc 2 souffre d'un Statement 
du point de conversion qui depend de Terreur sur y ef f, rnais cet etalement 
intervient toujours en sens oppose dans le cas ou Ton inverse un couple 
10 emetteur/recepteur donne. 

La coherence laterale, c'est a dire le deport AX entre image avant et 
arriere est done nulle lorsque la valeur de y e ff est correcte. 

Ainsi, la condition AX = 0 est equivalente a une determination 
exacte de y e ff ou a une determination exacte du deport Xc du point de 
15 conversion (determination dite de la configuration a point de conversion 
commun ou PCC). 

C'est ce qui est illustre sur les figures 3a et 3b. La figure 3a 
represente Tune en dessous de I'autre les images sismiques migrees avant 
et arriere obtenues pour y eff quelconque. Limage avant ne correspond a 
20 Hmage arriere qu'avec un decalage AX. La figure 3b represente quant a elle 
Tune en dessous de I'autre les images sismiques avant et arriere obtenues 
apres determination de la valeur de Yeff, c'est a dire apres mise en 
configuration a point de conversion commun. II n f y a alors pas de decalage 
entre les images avant et arriere. 
25 Le procede met done en oeuvre un traitement permettant de 

determiner la valeur de yeff pour laquelle le decalage AX est nul. 

Cette determination de la valeur de y e ff pour laquelle le decalage AX 
est nul se fait par exemple en posant 
vp Q = a vp 0 et vs p = p vs 0l 
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ou vpo et vso sont les valeurs approximatives de vitesses reelles locales de 
compression et de cisaillement determinees au moyen de migrations mises 
en ceuvre sur les acquisitions PP et PS initiates 
et ou a et p sont des variables. 

On peut par exemple utiliser pour vp 0 et vs 0 des valeurs globales 
correspondant a une estimation d'une vitesse moyenne a travers toutes les 
couches, y compris la couche analysee. On peut montrer qu'il est possible 
d'utiliser, au lieu des invariants classiques Vp, Vs, Tp, Ts, Vc, y e tf et y 0> les 
invariants Tp, Fp, Ts, Fs, ou encore Tp, Fp, Tc, Fc t ou Fp, Fs et Fc v6rifient 

Fp = Tp.Vp 2 

Fs = Ts.Vs 2 

Fc = Tc.Vc 2 

Ces invariants sont en effet mathematiquement equivalents aux 
invariants classiques Vp, Vs, Tp, Ts, Vc, yeff et y 0 . lis ont Pavantage d'etre 
facilement calculable a partir d'un champ de valeurs vp, vs, puisque que Ton 

a : 



ety e ff=Fp/Fs 

ou vp et vs sont les valeurs de vitesses reelles locales de 
compression et de cisaillement, ou I designe la coordonnee de la 
profondeur dans le sous-sol, ou z est la valeur de cette coordonnee de 
profondeur a la surface inferieure de la derniere couche d analyser et ou zo 
est la valeur de cette coordonnee de profondeur a la surface superieure de 
cette couche ou d'une couche au dessus de celle-ci, 

Par ailleurs : 

Tp, Fp, Ts, Fs constituent les quatre parametres independants 
fondamentaux dont tous les autres se deduisent. 
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Yeff = Fp/Fs, 
yo = Ts/Tp, 
Fc = Fp + Fs, 
Tc = Tp + Ts, 

5 Vp= yl(Fp/Tp) 

Vs = yj(Fs/Ts) 

Vc= yj(Fc/Tc) 
etc. 

II en resulte que les combinaisons quatre par quatre suivantes sont 
1 0 suffisantes pour determiner completement le probleme. 
(Tp, Fp f Ts, Fs) 
(Tp, Fp. Tc, Fc) 

(FC, TC, yo Yeff) 

(ce sont les plus utiles), 
15 et encore : 

(Tp, Vp, Ts, Vs) 
(Tp, Vp, Tc, Vc) 
etc. 

Pour chaque couple (a, p), on remplace les parametres 70. 7eff. Tp, 
20 Fp, Fc par 

yo' = o/p * 70 
7eff = a/p * 7eff 
Tp' = TpM1+Yo)/(1+Yo') 
Fp' = Fp* a 2 *(1+Yo)/(1+Yo') 
25 Fc 1 = Fc * ap * (1+y 0 )/(1+ Yo') * (1+ Yeff')/(1 + Yeff) 

et on met en oeuvre sur les donnees sismiques une migration PS 
correspondant a ces nouveaux parametres. 

Apr6s migration des donn6es sismiques, on calcule la correlation 
entre les images sismiques avant et les images sismiques arriere. 
30 On choisit le couple (a, p) pour lequel cette correlation est la plus 

importante. 
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Ce traitement est avantageusement mis en oeuvre en utilisant p = 
1/<x, ce qui permet de n'avoir d faire varier que la variable a (contre 
factorisation de vs et vp). 

Une fois y e ff ainsi determine, on met en oeuvre dans une nouvelle 
5 etape (etape 3 - focalisation H, la lettre H d6signant en effet classiquement 
le deport source-recepteur) un traitement destine £ permettre de raffiner la 
valeur du parametre Fc. 

A cet effet. on fait varier vp et vs tout en maintenant y e if constant. 
Par exemple, on pose vp a = a vpi et vs a = a vsi , oil vpi et vsi sont 
10 les valeurs determinees pour vp et vs dans T6tape 2, et ou a est une 
variable que Ton fait varier (co factorisation de Vp et Vs). 

On determine les nouveaux parametres y 0 \ Yeff 1 , Tp\ Fp\ Fc f 
correspondant a ces vitesses vp a et vs a (en rempla$ant p par a dans la 
formulation donnee ci-dessus pour ces parametres) et on met en oeuvre un 
15 traitement de migration en utilisant une moderation qui correspond a ces 
parametres. La migration produit cette fois-ci des collections a point de 
conversion (candidat) commun. Au lieu de regarder une image pour un 
deport avant et une image pour un deport arriere sur toute une ligne image, 
on regarde cette fois-ci en un ensemble de positions spatiales discretes 
20 (position PCC ou Point de Conversion Commun) alignees selon la direction 
de d6port consid6ree. 

La valeur retenue pour Fc est celle qui correspond a la valeur de a 
pour laquelle Talignement selon la direction de deport est optimale. 
L'analyse se fait sur autant de positions CCP que souhaite. 
25 yeff et Fc ayant ainsi ete determines a Tissue de cette troisieme 

6tape, il en va de meme pour Fs et Fp. 

II reste neanmoins a determiner Tp et yo, qui dependent Tun de 

I'autre. 

y 0 commande la condition de mise en profondeur commune des 
30 images PP (mais ne suffit pas d I'assurer sauf dans Thypothese incertaine 
d'une parfaite isotropie pour les ondes P et S). A contrario, une mesure 
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ind6pendante de yo. contrainte par la condition de mise en profondeur 
commune PP et PS donne un acces d Panisotropie du milieu. 

Tp est avantageusement determine, dans une etape 4, a partir du 
champ de vitesse vp determine a partir de I'analyse des ondes PP. On 
5 notera que la determination de Tp n'est pas indispensable a la 
determination des parametres Fp et Fs qui ont 6te determines de fa?on 
coherente grace au traitement de l'6tape 2 sur les ondes PS. 

On dispose ainsi A Tissue de cette etape 4 de tous les parametres 
de modelisation PS et PP : Tp, Fp, Ts, et Fs sont connus. 
10 Dans une etape 5, on met en oeuvre un traitement de mise en 

profondeur commune des modeles de vitesses S et P. En effet, les 
determinations des parametres Tp, Fp, Ts et Fs faites dans les etapes 
prec6dentes ont 6te faites sur la derniere couche et ne prejugeaient pas de 
1 I'existence d'une eventuelle anisotropie. Toutefois, la determination de ces 
15 parametres invariants a une profondeur commune aux modeles P et S 
n£cessite de prendre en compte Tanisotropie (en effet les invariants T et F 
font intervener des vitesses locales differentes en cas d'anisotropie). 

On utilise a cet effet un traitement se fondant sur un modelisation 
de Tisotropie transverse verticale ("VTI" ou "Vertical axis Transverse 
20 Isotropy" selon la terminologie anglo-saxonne utilisee par I'homme du 
m6tier) du type de celle proposee par Thomsen dans : 

[3] - Weak elastic anisotropy ; Geophysics, 51, 1954-1966 - 
Thomsen, 1986. 

et utilisant notamment les parametres d'anisotropie 5 et a introduits 
25 par Thomsen dans sa publication pr6citee. N 

(On rappelle que ces parametres d'anisotropie verifient en 
particulier : 

V nmop = V n sJ\ + 2S et V nmm = V xQ yfulZ 

oil Vnmop et V nrnos sont les vitesses NMO apparentes de 
30 compression et de cisaillement et ou Vpo et Vso sont les vitesses verticales). 

La vraie profondeur Zr est commune aux vrais modeles de vitesses 
P et S, et est reliee a Fp, Tp dans le premier modele , et a Fs, Ts dans le 
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second, par les deux parametres d'anisotropie 5 et o independants entre 
eux. 

Or, ces parametres 5 et <y sont contraints, en termes de valeurs 
moyennes apparentes, par la relation 
5 Yefl ^o = (1+28)/(l+2a). 

De plus , les parametres invariants Tp, Fp, Ts et Fs etant des 
integrates sur la profondeur, ils se pretent naturellement a une analyse 
couche par couche ("effeuillage" ou « layer stripping »), le probleme se 
ramenant £ une determination locale des invariants au toit de la derniere 
10 couche. 

On va alors determiner 5 inconnues locales, (v Pf v Si Z, 5, a) - c'est a 
dire les vitesses correspondant e la propagation en mode P et a la 
propagation en mode S, la profondeur commune de la base de la couche, et 
les deux parametres d'anisotropie - avec quatre mesures, Tp, Ts, Fp, et Fs. 
15 Une determination complete requiert soit de Tinformation a grand deport, 
soit de llnformation de puits pour contraindre Z ou 5. 
Par exemple, on impose 8 = 0. 

La focalisation et le calage des images PP et PS revient alors a 
resoudre le systeme d'equations suivant pour en trouver les solutions Z, v p , 
20 v s , a 

Fs = £vs.J\ + 2crxil 

25 

Dans une derniere etape (6tape 6), on met en ceuvre un traitement 
de courbure a grand deport. La resolution de I'analyse de vitesse des ondes 
P et S en milieu anisotrope, par des approximations isotropes et court 
deport n'a plus laisse qu'un seul degre de liberte au probleme: le 
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comportement a grand d6port ou I'anisotropie se manifeste le plus 
clairement. Le comportement des temps de trajet d grand deport est sonde 
par migration avant addition en faisant varier les parametres d'anisotropie 
tout en maintenant constant le rapport (1+28)/(1+2o). 
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REVEND1CATIONS 



1. Precede de prospection sismique selon lequel on emet dans le 
sous-sol une onde sismique de compression et on recueille a Paide de 
5 capteurs des donnees sismiques ayant au moins une composante de 
cisaillement et selon lequel on traite les donnees correspondant a cette 
composante de cisaillement pour en deduire une information sur la geologie 
du sous-sol, caracterise en ce qu'on determine une estimation du rapport 

£ vp.dl y/ £ vs.dl , ou vp et vs sont les valeurs de vitesses r6elles locales 

10 de compression et de cisaillement, ou I d6signe la coordonnee de la 
profondeur dans le sous-sol, ou z est la valeur de cette coordonnee de 
profondeur a la surface inferieure de la derniere couche a analyser et ou z 0 
est la valeur de cette coordonnee de profondeur a la surface superieure de 
cette couche ou d'une couche au dessus de celle-ci, 

15 et on inverse les donnees sismiques afin de deduire les valeurs de vitesses 
locales de compression et de cisaillement pour ladite couche a analyser, en 
utilisant une moderation dans laquelle cette estimation est utilisee pour le 
parametre invariant y e ff. 

2. Procede selon la revendication 1, caracterise en ce qu'on 

20 determine le paramdtre y e ff en mettant en ceuvre, pour differentes valeurs 
possibles pour celui-ci, un traitement de migration des donnees sismiques 
qui correspondent a la composante de cisaillement, et en d6terminant la 
valeur du parametre y e ff pour laquelle les images sismiques avant et arriere 
sont les mieux correl6es. 

25 3. Proc6de selon la revendication 2, caracterise en ce que pour 

faire varier le parametre yeff, on pose 
vp Q = a vpo et vs p = p vs 0l 

ou vpo et vs 0 sont les valeurs approximatives prealablement 
determinees pour vp et vs, 
30 et on fait varier a la fois les variables a et p. 
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4. Procede selon Tune des revendications 1 £ 3, caracterise en ce 
que la modelisation fait intervenir comme parametres invariants au moins 
quatre des parametres y 0 , y e ff. Tp, Fp, Tc, Fc ? avec y 0 = Ts/Tp, y eff = Fp/Fs, 
Tc = Tp + Ts oil Tp et Ts represented respectivement des temps de trajets 
5 verticaux des ondes de compression et de cisaillement, ou Fp est tel que 
(Fp/Tp) 1/2 represente une vitesse de compression et ou Fc est tel que ((Fc- 
Fp)/Ts) 1/2 represente une vitesse de cisaillement. 

5. Procede selon les revendications 3 et 4 prises en combinaison, 
caracterise en ce que, lorsque Ton fait varier les variables a et p, on 

10 remplace les parametres y 0 , y e ff. Tp, Fp, Fc par 
yo' = o/p*y 0 
Yeff' = Ct/P * yeff 

Tp' = Tp*(1+Yo)/(1+yo') 
Fp' = Fp* a 2 *(1 + Yo)/(1+yo') 
15 Fc l = Fc * ap * (1 +y 0 )/(1 + yo f ) * (1 + yeff)/(1 + yeff) 

et on met en oeuvre sur les donn6es sismiques une migration 
correspondent a ces nouveaux parametres. 

6, Procede selon la revendication 5, caracterise en ce que pour faire 
varier le parametre y eff , on pose p = 1/a et on fait varier a. 

20 7. Procede selon Tune des revendications precedentes, caracterise 

en ce qu f a Tissue de la determination du parametre y e ff, on fait varier vp et 
vs tout en maintenant y e fr constant et on determine le parametre Fc pour 
lequel Talignement selon la direction de deport est maximal. 

8. Procede selon la revendication 7, caracterise en ce pour faire 

25 varier vp et vs, on pose vp a = a vpi et vs a = a vs 1t ou vpt et vs! sont des 
valeurs determinees pour vp et vs dans Tetape 2, et on fait varier la variable 
a. 

9. Procede selon Tune des revendications 7 ou 8, caracterise en ce 
qu'3 Tissue de la determination du parametre Fc, on determine le parametre 
30 Tp et/ou le parametre y 0 = Ts/Tp. 
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10. Proc6de selon la revendication 9, caracterise en ce que le 
param6tre Tp est avantageusement d6termin6 & partir du champ de vitesse 
vp d6termin§ a partir de t'analyse de la composante de compression des 
donnees sismiques. 

5 11. Proc6d6 selon Tune des revendications 9 et 10, caracterise en ce 

qu'on met ensuite en oeuvre un traitement de mise en profondeur commune 
des modeles de vitesses S et P. 

12. Procede selon la revendication 11, caracterise en ce qu ! a Tissue 
du traitement de mise en profondeur, on met en ceuvre un traitement a 
10 courbure a grand d6port en faisant varier les parametres d'anisotropie 6 et a 
tout en maintenant constant le rapport (1+28)/(1+2o). 
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